Superfluidity of the BEC at finite temperature. 
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We use the classical fields approximation to study a translational flow of the condensate with 
respect to the thermal cloud in a weakly interacting Bose gas. We study both, subcritical and 
supercritical relative velocity cases and analyze in detail a state of stationary flow which is reached 
in the dynamics. This state corresponds to the thermal equilibrium, which is characterized by the 
relative velocity of the condensate and the thermal cloud. The superfluidity manifests itself in the 
existence of many thermal equilibria varying in (the value of this velocity) the relative velocity 
between the condensate and the thermal cloud. We pay a particular attention to excitation spectra 
in a phonon as well as in a particle regime. Finally, we introduce a measure of the amount of the 
superfluid fraction in a weakly interacting Bose gas, allowing for the precise distinction between the 
superfluid and the condensed fractions in a single and consistent framework. 

PACS numbers: 03.75.Hh, 03.75.Kk, 47.37.+q 



I. INTRODUCTION 

The tendency of liquid helium-4 cooled below 2.19 K to 
flow without any apparent friction is known for decades 
■ This superfluid behavior is commonly associ- 
ated with the phenomenon of a Bose-Einstein condensa- 
tion 0. The connection has been mostly based on the 
similarities in the nature of both phenomena, as both 
are a large-scale manifestations of the quantum nature. 
However, the detailed theoretical investigation of super- 
fluidity and condensation in liquid helium 5] is a great 
challenge, as it is a strongly interacting quantum system 
and cannot be effectively described with a mean-field or 
perturbative approaches. The present-date understand- 
ing of superfluidity comes from the paper of L. Landau 



Experimental achievement of Bose-Einstein conden- 
sates in dilute trapped atomic gases Q, 0, gives a 
unique possibility to study their superfluid properties un- 
der condition of very weak interactions, ft has opened 
an opportunity for revisiting the concept of superfluid- 
ity. In this context the issue of the atomic Bose-Einstein 
condensate under rotation has attracted great interest 
[Ifj| . It has been demonstrated experimentally that Bose- 
Einstein condensate circulates around quantized vortex 
lines. In experiments with gasous BEC the existence 
of these vortices and the properties of scissor modes 
[TTl Il2l IT3L IT3. PHI are considered the main manifes- 
tations of superfluidity. On the other hand, a little is 
known about superflui dity in a translational motion of 
the atomic condensate [l6|, [Tt], EH • There are still no ex- 
periments that would correspond to the demonstration of 
frictionless non-rotary flow, which is the most intuitive 
manifestation of this phenomenon. 

Dilute atomic Bose gases are simpler to describe than 
liquid helium, with many succcsful implementations of 
mean-field theories. It is a purpose of this paper to 
present a simple model of superfluidity in a weakly inter- 



acting Bose gas using the classical fields approximation. 
The classical fields approximation is the non-perturbative 
mean-field model of a Bosc-Einstcin condensate at fi- 
nite temperatures. It turned out to be very success- 
ful in describing dynamical and thermodynamic prop- 
erties of atomic condensates and their excitation spectra 

d US HIE! IHH]. :t all0WS a descri P tion of the con- 
densate and a thermal cloud treating both components 

on equal footing without any arbitrary splitting of the 

system into two parts. 

In Section |n] we briefly describe the classical fields ap- 
proximation. For ease of calculations we restrict our- 
selves to a simplified geometry of 3D box with peri- 
odic boundary conditions. Then we prepare initial states 
which will allow us to study the relative motion of the 
condensate and the thermal cloud. 

In Section IlIII we present numerical results. We find 
that the flow is indeed frictionless below the Landau crit- 
ical velocity c and that the stationary state depends on 
yet another control parameter - namely the relative ve- 
locity between the condensate and the thermal cloud. 
We believe that in this paper we are extracting for the 
first time the basic properties of the superfluid flow in 
direct numerical simulation of the dynamics. We ob- 
tain the excitation spectra of such a superfluid system 
and observe the equipartition of energies in quasiparti- 
cle modes, which has been previously found in case of 
a stationary BEC [l^. 124^ . We also notice pulling of 
Bogoliubov quasiparticles by the condensate during the 
thermalization period. 

In Section ll VI we analyze the properties of Bogoliubov 
quasiparticles. We derive analitically excitation spectra 
of the system, obtaining results that agree with our nu- 
merical simulations. We analyze the transformation of 
quasiparticle momentum under Galilean transformations 
and find that it behaves as a particle of mass m, simi- 
larly to atoms in the system and unlike acoustic waves. 
This allows us to directly relate superfluid and condensed 
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fractions, what produces surprisingly simple formula in 
the case of 3D box. We base our reasoning on the ob- 
servation that Bogoliubov quasiparticles are the normal 
modes of the system and we associate with them the 
normal (non-superfmid) fraction, similarly to Landau. In 
Section El we conclude. 



II. THE METHOD 

We consider a weakly-interacting Bose gas consisting 
of N atoms confined in a 3D box potential with periodic 
boundary conditions. The atoms interact via a contact 
potential V(r — r') — 47r ^ a S^(r — r'), where a is the 
s-wave scattering length. They can be described by the 
second-quantized Hamiltonian: 
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(2.1) 

where L is the size of the box, m is the mass of an 
atom and \1/ is a bosonic field operator satisfying equal 
time bosonic commutation relation [^f{r,t),"^f'(r',t)] — 
S^(r — r'). The corresponding Heisenberg equation for 
^> is of the form: 
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We expand the field operator in the basis of plane waves: 
^ = -^== ^ k e 27 ™ kr / L dk(i), where annihilation operators 
Sk destroy particle in mode k and satisfy a commutation 
relation [ak, a^,] = <5k,k'- 

Now we apply the classical fields approximation |2lLl23 . 
12311251 l26j . For all modes that are occupied by sufficiently 
large number of atoms we subtitute annihilation opera- 
tors with c-numbers: &k — > s/Na^. After neglecting all 
remaining operator terms we obtain a set of equations: 
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We choose wavectors of "classical" modes ak to span 
a grid of size 32x32x32. They compose the normalized 



mean-field wave function ip 
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This wavefunction obeys the Gross-Pitaevskii equation 
on a spatial grid: 
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where g — aN/itL is the interaction strength and e = 
4ir 2 h 2 /mL 2 is our unit of energy (the corresponding unit 
of time is h/e). 

The wavefunction if) represents here both the conden- 
sate and the thermal cloud, contrary to the typical in- 
terpretation of the GP equation where ip is just the con- 
densate wavefunction. A measurement, which introduces 
a coarse-graining is essential for distinguishing between 



these two fractions. Populations of various modes can 
be extracted from diagonalization of the time-averaged 
density matrix j2^, H3| and the dominant eigenvalue rep- 
resents the condensed fraction. The time- aver aging de- 
stroys the coherence between the modes and leads to a 
mixed state out of a pure state \tp >< The wave- 
function i[>(t) can be interpreted as a single experimental 
realization of evolving system. 

To obtain a stationary solution corresponding to non- 
zero temperature we start with the ground state wave- 
function and randomize its phase. The system evolves 
towards a thermal equilibrium, which depends only on 
the total energy per particle i?, the interaction strength g 
and the length of the box L |23| ■ The condensate appears 
in the zero momentum mode and atoms in the remaining 
modes compose a thermal cloud. The dependence of the 
resulting state on the interaction strength g = aN/irL 
instead of a and N separately results in an ambiguity. In 
order to assign the temperature T to our system in a way 
independent of the choice of the numerical grid we use a 
scheme developed in [2l| ; we assign temperature to tem- 
perature of the ideal gas that has the same condensate 
fraction. 

In the considered case of the box potential eigenmodes 
of the stationary single-particle density matrix can be 
only composed of modes with the same |k|, due to sym- 
metry constrains. Thus the plane waves set a natural 
basis. As it turns out, modes k combine with modes 
-k to create Bogoliubov quasiparticles of amplitude <5k, 
which evolve with a single dominant frequency (the rala- 
tion between and ak will be derived in Section ITV|l . 
Bogoliubov quasiparticles are the normal modes of the 
system |19|: for high momenta, when the kinetic energy 
is much larger than the interaction energy, they become 
indistinguishable from atomic modes ak- 
in a stationary case the central energies of quasiparti- 
cles are given by [13: 



£k = A* + e\/w(k) 2 - (gn c ) 2 



(2.5) 



where w(k) = + gn c is introduced for simplicity and 
n c is the condensate fraction. \x is the frequency of the 
condensate mode and has an interpretation of the chem- 
ical potential. Similarly to the case of He-II they have 
a phonon-like spectrum at low momenta (see Fig. 5 in 
[ljj). In the classical fields approximation quasiparticles 
obey the equipartition of energies [n|> suggesting that 
the system is in a thermal equilibrium and not just in 
a steady state. Note, however, that the condensate it- 
self is not populated according to the equipartition. For 
parameters considered in our paper the Bose statistics 
leads to the equipartition for the considered Bogoliubov 
quasiparticles, suggesting that thermodynamics should 
be studied in quasiparticle and not atomic modes. 

The solutions discussed here represent the situation in 
which the condensate does not move with respect to the 
thermal cloud. In order to study the possible superfluid 
flow we need to put the condensate in motion. In our nu- 
merical simulations we take a stationary state described 
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above and apply a momentum kick to the condensate. 
The remaining atoms are pushed in the opposite direc- 
tion, so that the total momentum of the system is close 
to zero We can also apply the momentum kick to 

any arbitrary number of modes, in particular to the con- 
desnate along with surrounding modes which merge into 
quasiparticles (phonon modes). 

In our simulations we take stationary solutions for ca. 
60000 atoms and the scattering length of 18 nm, confined 
in a box of lenght 5 /im, with the condensate popula- 
tions of 65% and 8%. The temperatures are 580 nK and 
1 fj,K, respectively. We apply the procedure described 
above, controlling the initial relative velocity v of the two 
fractions, by placing the condensate in a chosen mode 
k c . The resulting states are composed of the condensate 
(possibly with some surrounding atomic modes) and the 
remaining thermal cloud flowing in opposite directions. 
These states are arbitrarily prepared and no longer rep- 
resent an equilibrium of the system - they should experi- 
ence some thermalization. 

Can we observe any frictionless (superfluid) flow be- 
tween the two fractions? 



III. NUMERICAL RESULTS 

The steady state depends on the initial relative ve- 
locity v between the condensate and the thermal cloud. 
There exists a critical velocity c = 2^ y/gn c which sepa- 
rates two regimes; n c is the new equilibrium condensate 
fraction. For a system with 65% (8%) of the condensate 
initially the resulting n c is approximately 50% (1.5%) 
for velocities lower than c. For relative velocities v be- 
low the critical velocity c the condensate remains flowing 
with respect to the thermal cloud. No drag due to ther- 
mal atoms is experienced and the flow is superfluid (Fig. 

Note the initial thermalization. We have observed 
that during this transient period the phonon modes are 
being dragged along with the condensate. Initially push- 
ing low-excited modes along with the condensate reduces 
this time and increases the equilibrium condensate frac- 
tion (see inset in Fig. [JJ. 

Above the critical velocity the condensate experiences 
a drag from thermal atoms and slows down. The conden- 
sate population slips to a slower moving mode (Fig. 0). 
Close to the critical velocity the system is unstable. Slips 
will occur at different times, depending on the particular 
realization of the initial conditions (Fig. [3J). Bundling 
phonon modes with the condensate results in stabilizing 
the near-critical flow. This is due to an increase in the 
resulting equilibrium condensate fraction which shifts the 
critical velocity slightly upwards. 

The excitation spectrum of quasiparticles of the sys- 
tem with the condensate in k c mode has been plotted in 
Figure 01 It can be well described by (see below): 

e k = ft + evMk - k c ) 2 - (gn c ) 2 + e(k - k c )k c . (3.1) 
Relative energies of quasiparticles moving faster than the 
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FIG. 1: Fraction of the condensate versus time for the case 
of superfluid flow with initial velocities between the thermal 
cloud and the condensate equal v = 0.48c and v — 0.82c. 
The equilibrium condensate fraction is 51.6% and 49.6%, re- 
spectively. The initial states were prepared by pushing the 
condensate without phonon modes. In this case dragging 
of phonons has occured. The inset shows a situation where 
phonon modes have been initially pushed with the condensate 
and v = 0.79c. The resulting condensate fraction is 53.4% in 
this case. 

condensate are shifted upwards, and moving slower than 
the condensate are shifted downwards with respect to the 
stationary state, due to the Doppler-like shift (the last 
term in Formula 13. 1|) . Note, that numerically obtained 
spectra of quasiparticles have finite widthjwhich can be 
interpreted as the their finite life-time [jjj- Still, the 
quasiparticles obey the equipartition of energies. Their 
populations ris k = <5£<5k and energies ek fulfill a formula: 

ns k {ek - (Jt) = k B T, (3.2) 

as can be seen in Figure [5^. This indicates that the 
system undergoing a superfluid flow is in a thermal equi- 
librium. However, the equipartition occurs not in the 
condensate nor the thermal cloud's frame, but in the 
center-of-mass frame instead. 

Please confront the relation 13 . 21 with a similar one ob- 
tained for atomic populations = aj^evk (Fig. Et>), 
which do not obey the equipartition. For phonon modes, 
where quasiparticles differ significantly from atoms, their 
populations are smaller than relative atomic populations. 

For starting n c = 8% the critical velocity corresponds 
to a momentum lower than the momentum of the k = 
(1, 0, 0) mode and thus any motion of the condensate with 
respect to the thermal cloud results in a friction. 

Within our accuracy the observed critical velocity is 
identical with the one obtained from Landau's criterion 
applied to the excitation spectrum in the stationary 
case (Formula^)). Note, however, that even for subcrit- 
ical velocities there exist thermal atoms that move super- 
critically with respect to the condensate. Applying the 
Landau criterion to each quasiparticle separately would 
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FIG. 2: Population of the condensate versus time for a super- 
critical flow. The starting relative velocity of the condensate 
and the thermal cloud is v = 1.17c (a) and v — 1.75c (b). The 
mode initially occupied by the condensate becomes depleted 
due to friction with the thermal cloud (I) and the condensate 
reappears in a slower moving mode (II). Note, that in (a) the 
condensate is directly transferred to the neighbouring momen- 
tum state, where it flows with a stationary relative velocity 
v = 0.46c with respect to the thermal cloud. On the contrary 
in (b), where the initial velocity is very high, the flow is much 
more complicated. The condensate population smears into 
multiple modes to finally agregate in a mode whose relative 
velocity equals v — 0.52c. In the transient time many suc- 
ceeding slips may occur, as the condensate slows down gradu- 
ally. As a result of these slips the stationary relative velocity 
of the condensate and the thermal cloud is subcritical. The 
equilibrium condensate population is 47.9% in (a) and 42.9% 
in (b). 



always result in a net drag force, as there is more super- 
critical atoms moving opposite to the condensate than 
supercritical atoms moving along with the condensate. 

The manifestation of superfluidity in our model is the 
existence of many thermal equilibria, all varying in the 
relative velocity between the condensate and the thermal 
cloud. Such a velocity is an additional control parameter. 
Please note that we have not yet presented any quanti- 
tative measure of the amount of the superfluid fraction. 
In the next section we present analytical discussion and 
derive the formula for the superfluid population. 




FIG. 3: Fraction of the condensate versus time for a flow 
with the critical velocity. (I) The condensate has been pushed 
without phonon modes and the initial relative velocity be- 
tween the condensate and the thermal cloud is v = c. The 
two plots represent different realizations of the same setup. 
Slips occur at different time and show that the system is un- 
stable in the critical region. (II) Phonon modes have been 
initially pushed with the condensate and stabilize the flow, 
which is now superfluid. The relative velocity between the 
condensate and the thermal cloud is now v = 0.98 c due to a 
slightly greater condensate population. 




FIG. 4: The excitation spectra of quasiparticles in a direction 
along the momentum kick. Numerical results (stars) were 
obtained for the case of the initial condensate fraction n c = 
0.65 pushed together with the phonon modes with the initial 
velocity relative to the thermal cloud equal v = 0.82c. Dots 
depict a fit with Formula 13.11 The equilibrium condensate 
fraction is n c = 0.534. Spectra taken in directions different 
than presented experience smaller tilting in agreement with 
Formula EH1 



IV. ANALITICAL RESULTS 

We consider here Bogoliubov excitations at non-zero 
temperature. The starting point of our investigation will 
be Eq. 12.21 for the Bose field operator ^S>. Let us choose 
a reference frame in which the condesate is in k c mode. 
Changing k c will be interpreted for a while as changing 
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the reference frame. We decompose if? into a basis of 
plane waves: 



(4.1) 



with indices chosen in such a way that a' represents the 
anihilation operator of the condensate mode. The equa- 
tion for a' k is: 



An ha 
mL 3 



E 

ki, k 2 



a U 1 a k 2 a k+k 1 -k 2 



(4.2) 



It is a non-linear operator equation. We linearize it in 
a way similar to the Bogoliubov-Popov approximation. 
For certain elements in the sum in Eq. 14.21 we group 
two operators in pairs and replace them with apporpri- 
ate mean values. And thus, terms a^ak are replaced with 
< a k <2k >= Nrik, which do not depend on time for a sta- 
tionary state. a' a' is also replaced by iVn c e~ 2lAlt , where 
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FIG. 5: Energies localized in quasiparticle modes (a) and 
atomic modes (b) versus their momenta for parameters as 
in Fig. [I] Squares in (b) are analitical predictions for ener- 
gies in atomic modes based on Formula 14. 121 while diamonds 
are numerical results. The difference between analitical and 
numerical results comes from neglecting the finite width of 
the quasiparticle spectra in our analitical calculations. 



the chemical potential /i is obtained from the equation 
for a' Q . The remaining terms in the summation, which 
are still non-linear in operators, are neglected. Results 
from ^3 suggest that this approximates the full dynam- 
ics rather well. The neglected terms become important at 
temperatures close to the critical temperature, limiting 
the scope of our analysis. 

The approximate equation for a' is: 



ihd t a' (t) = e-± a' + e(gn c + 2gn T )a' , (4.3) 

where n-r is the thermal fraction, defined as nx = 1 — 
n c . It gives the leading contribution to the frequency 

k 2 

of the condensate mode u = gn c + 2griT + -f , which is 
in agreement with results from the Popov approximation 
|28j . Corrections to a have been discussed in 0. For 
k 7^ we obtain: 



ihdta' k (t) = e 
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where: 



e+ = ey/iv jk) 2 - (gn c ) 2 + ekk c 
e k = -e v /w(k) 2 - (gn c ) 2 + ekk c . 



(4.5) 



(4.6) 



Now we perform Bogoliubov transformation. We in- 
troduce new bosonic anihilation operators: 



lf k a k + ViaV~ 2l(AI+£kkc)t - 



(4.7) 



which destroy a Bogoliubov excitation in mode k. The 
commutation relation [Sk, = <$k,k' implies that |f| 2 — 
\V\ 2 = 1. We look for normal modes of the system (which 
oscillate with a single frequency only), which are linear 
combinations of solutions 14.51 The resulting Bogoliubov 
quasiparticles evolve with the frequency fi + e k : 

S k (t) = i k (0)e-'("+* t (4.8) 
and the coefficients Uk and Vk are given by: 



u2 _ (^(k) + ^(k)2-(g, l<: )2) 2 



(4.9) 



(w(k)+V"(k) a -(flnc) 2 ) -(S«c) 2 
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Energies form the excitation spectrum of the system. 
They agree with values obtained from the numerical sim- 
ulations. However the calculated value of the chemical 
potential \x agrees rather with the Formula 14 of 0] 
than with the simplified Equation 14.31 The difference, 
however, shifts the whole energy spectrum only. Note, 
that this simplified approach ignores the finite width 
(which determines the life-time) of the quasiparticle spec- 
tra, which we observe in our numerical simulations. 

We can also perform an inverse Bogoliubov transofma- 
tion. It is: 

a k = C/ k <S k - ViiL k e- 2 ^ +ekk ^ (4.10) 
and the resulting atomic populations are: 

Nn k =< a£a' k >= u£ < 5£<S k > +V k < 5_ k 5l k > . 

(4.11) 

Note that U k = f/_ k , V k = F_ k . 

The resulting Bogoliubov quasiparticles <5 k have mo- 
menta ps k — h(k + k c ) and their coefficients C/ k and 
are independent of the choice of the reference frame (k c ). 
Thus they transform under Galilean transform as parti- 
cles of mass m, similarly to atoms in our system. 

Let us now switch to the semi-classical model again, 
by substituting c-numbers for the operators (a' k — * a k 
and (5 k — * <5 k ). By doing this we neglect the antinormal 
ordering of <5 k in Formula 14.111 which accounts for the 
quantum depletion. On the other hand we can relate 
the results to our numerical simulations, allowing us to 
calculate the amount of superfluid fraction. 

In this semi-classical approach Formula 14. 1 II takes the 
form: 

n k = Uln 5 ^ + t£n*_ k , (4.12) 

where ns k = ^ k <5 k are quasiparticle populations. Assum- 
ing we know quasiparticle populations we can calculate 
populations of atomic modes with this formula. They 
are greater than the corresponding quasiparticle popula- 
tions due to U k > 1. These analitical calculations agree 
with our numerical results (see Fig. \Sjp for comparison of 
atomic populations). 

In order to identify the superfluid part we consider the 
total momentum P of the system. It can be naturally 
expressed as a sum over momenta of individual atoms: 

P = HnI n c k c + J2 «k(k + k c ) ] =P C + P T , (4.13) 

where atoms contribute either to the momentum of the 
condensate P c or to the momentum of the thermal cloud 
P T (for k^ 0). 

On the other hand quasiparticles have a well-defined 
momenta too and they are also the observed excitations 
that can be measured experimentally [2^, |3(], 0OJ . Their 
agregate momentum is Ps — NhJ2u=£o n $kO i + kc) and 
we can express the total momentum P: 

P = P Sk +P SF - (4.14) 



An additional momentum P$f appears here to account 
for the conservation of momentum. Because Ps and P 
are physically observable Psf should have some physical 
meaning. Applying the inverse Bogoliubov transforma- 
tion (|4.12|1 to atomic populations and using the relations 
for C/ k and one can show that: 

n T = n s + 2 ^ n Sk V^, ^ ^ 

This implies that: 

P SF = Nh(n c + n T - n s )k c , (4.16) 

where ns = X) k ^o ns k i s the total number of quasiparti- 
cles. Only in the condensate reference frame this mo- 
mentum is zero. It appears to come from a fraction 
nsF — n c + tit — ns — 1 — ns moving with the con- 
densate velocity. This mass represents condensed atoms 
along with a part of thermal cloud " freezed out" by the 
interaction. 

We have decomposed the system into two parts - one 
composed of excitations of the system and another one 
moving with the condensate speed. Applying the Landau 
criterion to such a system we conclude that the usf — 
1—ns is indeed the superfluid fraction, as indicated by the 
subscript we used, and Bogoliubov quasiparticles form 
the normal fraction in the system (ns). 

In other words the superfluid fraction is a sum of the 
condensate and atoms which do not take part in forming 
quasiparticles (which are thermally distributed). Thus 
atoms occupying highly excited modes do not contribute 
to the superfluid part, as for these modes quasiparticles 
are almost the same as atomic modes. Note that the su- 
perfluid fraction adds up with quasiparticles to the total 
number of atoms in the system, even though Bogoliubov 
quasiparticles do not obey the conservation of their num- 
ber. 

Due to the bosonic nature of Bogoliubov quasiparti- 
cles the superfluid fraction is greater than the condensate 
fraction and the difference increases with the interaction 
strength g. In our numerical simulations it is for exam- 
ple nsF = 60.4% compared to n c — 53.4% for relative 
velocity v = 0.82c between fractions. 

We can apply the same reasoning without restricting 
ourselves to the semi-classical description of atoms and 
quasiparticles. The same scheme applied to a system 
with quantum corrections taken into consideration shows 
that even at zero temperature the superfluid fraction is 
100% even though the condensed fraction is smaller than 
that. Because the quantum depletion originates from the 
reversed order of S k in the inverse Bogoliubov transform, 
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the formula for superfluid population takes the form: 
Nn SF = Nn c + J2 (&k < *£& > +V£ < <L u ^ k >) - 

E < > 

k^O 

= Nn c + Nn T - Nn 5 + N QDl 
(4.17) 

whereas the total number of atoms is given by: 

N = Nn c + J2 (Ul < Sl6 k > +V£ < <L k <f k >) 

k^O 

= Nn c + Nn T + N QD . 

(4.18) 

In these formulas Nng = X^k^o < ^k4 >> Nut = 
Nn$ + 2^ k ^ ^< i5 k (5k > V^j and the quantum deple- 
tion Nqo — Jik/o * s accounted separately. At zero 
temperature, when < <5 k <5k >= 0, both tit and ng vanish 
and the superfluid population equals the total number of 
atoms in the system. 

Even though such a basic model suffers from a diver- 
gence in quantum depletion terms resulting from the zero 
range potential, it clearly indicates that the quantum de- 
pletion belongs to the superfluid fraction. 

V. CONCLUSIONS 

In this paper we have presented a model of superflu- 
idity of a weakly-interacting Bose gas within the classi- 
cal fields approximation. In the direct numerical sim- 
ulation of the dynamics we have shown that below the 
Landau critical velocity c = f^jj-y/S^c the relative flow 
of the condensate with respect to the thermal cloud is 
superfluid and that the resulting state is in a thermal 
equilibrium. The relative velocity of this superflow is an 
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additional control parameter characterizing the thermal 
equilibrium. 

We have obtained the excitation spectrum of the sys- 
tem and we have analyzed in detail properties of Bo- 
goliubov quasiparticles, obtaining an agreement between 
numerical and analitical calculations. 

The superfluid fraction has been identified from the 
momentum transformation rules and an application of 
the Landau criterion. It turns out to be the conden- 
sate plus atoms from the excited modes which do not 
participate in the collective excitations of the system. 
These atoms can be viewed as bound to the condensate 
by the interactions. The results provide analitical for- 
mula for a difference between superfluid and condensed 
parts at any temperature and within a consistent and 
non-perturbative description. They also show that one 
does not require a strongly-interacting Helium-II, with 
its roton part of the spectrum and essential many-body 
interactions, to make this difference significant. 

In our approach it is also evident that quantum deple- 
tion occuring in a weakly-interacting Bose gas is also a 
part of the superfluid fraction, so that even at OK there is 
a difference between superfluid and condensed fractions 
- again similarly to He-II. 

The temperature dependent Bogoliubov-Popov spec- 
trum, together with our approach to superfluidity can 
form a basis for equilibrium statistical physics of inter- 
acting Bose gas at non-zero temperature. Such a study 
will be presented elsewhere [3^ |. 
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